1. Data Exploration

This week’s data analysis mostly deals with the human-dataset. The data contained a number of variables mapping human development and gender equality in 195 countries. This data was then cut down to eight variables. They are the following variables measuring development and gender equality in 195 countries:
+ Life.Expectancy refers to a citizens average life expectancy at birth
+ Education.Years is the amount of years a citizen is planned to spend in education
+ GNI.Per.Capita is the Gross National INcome adjusted for population
+ Maternal.Mortality is a ratio of deaths per 100,000 births
+ Teen.Birth is the number of births per 1,000 women ages 15 to 19.
+ Parliamentary.Participation is the ratio of women to men in parliament
+ Female.Secondary.Education is the percentage of females that attend secondary education
+ Female.Work is the percentage of females that participate in the labour force

Below are simple bar-plots graphically showcasing the variables used this week. As the graphs show, most of these variables do not follow a gaussian distribution, meaning that they are not bell-curves. The only exception to this is the variable mapping the planned years a citizen is supposed to spend in school.

Even further down, the reader can find the correlation plot, where "*" indicates statistical significance to the 0.05-level. As can be seen, three sets of variables emerge from the data:
+ Female.Work and Parliamentary.Participation, which are not strongly correlated with anything
+ Maternal.Mortality and Teen.Birth, which are strongly positively correlated only with one another.
+ The remaining variables on education, GNI, and life-expectancy, which are (semi) strongly positively correlated with each other and strongly negatively correlated with the preceding set of variables.

##  Female.Secondary.Education  Female.Work    Life.Expectancy Education.Years
##  Min.   :  0.90             Min.   :13.50   Min.   :49.00   Min.   : 5.40  
##  1st Qu.: 27.15             1st Qu.:44.45   1st Qu.:66.30   1st Qu.:11.25  
##  Median : 56.60             Median :53.50   Median :74.20   Median :13.50  
##  Mean   : 55.37             Mean   :52.52   Mean   :71.65   Mean   :13.18  
##  3rd Qu.: 85.15             3rd Qu.:61.90   3rd Qu.:77.25   3rd Qu.:15.20  
##  Max.   :100.00             Max.   :88.10   Max.   :83.50   Max.   :20.20  
##  GNI.per.Capita   Maternal.Mortality   Teen.Birth    
##  Min.   :   581   Min.   :   1.0     Min.   :  0.60  
##  1st Qu.:  4198   1st Qu.:  11.5     1st Qu.: 12.65  
##  Median : 12040   Median :  49.0     Median : 33.60  
##  Mean   : 17628   Mean   : 149.1     Mean   : 47.16  
##  3rd Qu.: 24512   3rd Qu.: 190.0     3rd Qu.: 71.95  
##  Max.   :123124   Max.   :1100.0     Max.   :204.80  
##  Parliamentary.Participation
##  Min.   : 0.00              
##  1st Qu.:12.40              
##  Median :19.30              
##  Mean   :20.91              
##  3rd Qu.:27.95              
##  Max.   :57.50

2. Pricipal Component Analysis on Non-Standardized Data

Below, Principal Component Analysis, relying on Singular Value Decomposition, is applied to the non-standardized “Human”-data explored above. This analysis is expected to turn out somewhat skewed, as PCA relies on variance to identify the Principal Components. Consequently, variables with larger average values will have an unreasonably high impact on the analysis.

## Standard deviations (1, .., p=8):
## [1] 18544.172057   186.283543    25.972416    20.074641    14.321634
## [6]    10.634338     3.721257     1.428190
## 
## Rotation (n x k) = (8 x 8):
##                                       PC1          PC2           PC3
## Female.Secondary.Education  -9.317458e-04  0.085169068 -0.3652346210
## Female.Work                  9.124960e-05 -0.031442122 -0.0180958993
## Life.Expectancy             -2.815825e-04  0.028224407 -0.0170864406
## Education.Years             -9.562916e-05  0.007556395 -0.0210421918
## GNI.per.Capita              -9.999828e-01 -0.005830843  0.0006412388
## Maternal.Mortality           5.655739e-03 -0.987500960 -0.1481355205
## Teen.Birth                   1.233962e-03 -0.125305410  0.9184673154
## Parliamentary.Participation -5.526462e-05  0.003211742 -0.0038388487
##                                       PC4           PC5           PC6
## Female.Secondary.Education  -0.8435797499 -3.770012e-01 -6.083913e-02
## Female.Work                 -0.3928157700  8.233860e-01  4.077784e-01
## Life.Expectancy             -0.0212996883  1.136966e-02 -6.669649e-02
## Education.Years             -0.0308785262  3.896982e-03 -2.437473e-02
## GNI.per.Capita               0.0002381021  6.651486e-06  2.062449e-05
## Maternal.Mortality          -0.0173448186 -3.955532e-02 -2.010447e-02
## Teen.Birth                  -0.3475453954 -1.381061e-01 -2.499459e-02
## Parliamentary.Participation -0.1075784134  3.989026e-01 -9.077136e-01
##                                       PC7           PC8
## Female.Secondary.Education   0.0303039622  3.118655e-02
## Female.Work                 -0.0093667016  6.654689e-03
## Life.Expectancy             -0.9901494274  1.161211e-01
## Education.Years             -0.1134676761 -9.925031e-01
## GNI.per.Capita               0.0001028639  1.871381e-05
## Maternal.Mortality          -0.0244043639 -7.101321e-04
## Teen.Birth                  -0.0127951595 -8.079546e-03
## Parliamentary.Participation  0.0704544742  1.925677e-02
## Importance of components:
##                              PC1      PC2   PC3   PC4   PC5   PC6   PC7   PC8
## Standard deviation     1.854e+04 186.2835 25.97 20.07 14.32 10.63 3.721 1.428
## Proportion of Variance 9.999e-01   0.0001  0.00  0.00  0.00  0.00 0.000 0.000
## Cumulative Proportion  9.999e-01   1.0000  1.00  1.00  1.00  1.00 1.000 1.000

As expected, the graph can tell us very little, and PC1 purports to explain (essentially) all of the variation. Let us standardize the data and try again.

3. Pricipal Component Analysis on Standardized Data

By standardizing the variables in the “Human”-data, the PCA is able to work on variables with comparable variances. Accordingly, variables on scales where high values are expected will not have a disproportionately high impact on the calculations. As can be seen below, the biplot produced by the standardized data can offer valuable insights into the data.

## Standard deviations (1, .., p=8):
## [1] 2.1359720 1.1170602 0.8767066 0.7202736 0.5530904 0.5299143 0.4493522
## [8] 0.3372776
## 
## Rotation (n x k) = (8 x 8):
##                                     PC1         PC2          PC3         PC4
## Female.Secondary.Education  -0.39874523 0.101737848 -0.210481695 -0.31033900
## Female.Work                  0.14439668 0.682158952 -0.575062704 -0.27300332
## Life.Expectancy             -0.43080653 0.003997077  0.073302641 -0.02783772
## Education.Years             -0.41858363 0.138149663 -0.073869337 -0.07719277
## GNI.per.Capita              -0.33914119 0.098797891 -0.338769060  0.84020987
## Maternal.Mortality           0.41991628 0.123208094 -0.145471957  0.28520785
## Teen.Birth                   0.40271826 0.088902996 -0.003213286  0.11830579
## Parliamentary.Participation -0.07626861 0.687286162  0.691544283  0.14537131
##                                     PC5         PC6         PC7         PC8
## Female.Secondary.Education   0.43437742 -0.50478735  0.48033587  0.12578655
## Female.Work                 -0.22413576  0.23657515  0.01076780 -0.04754207
## Life.Expectancy              0.04920168  0.57970641  0.04122211  0.68415054
## Education.Years              0.30831448 -0.11017948 -0.81713841 -0.13919229
## GNI.per.Capita              -0.01249472  0.01564322  0.16310711 -0.16583233
## Maternal.Mortality           0.05493343 -0.43780233 -0.24286926  0.67254029
## Teen.Birth                   0.81248359  0.36928738  0.07464930 -0.11761127
## Parliamentary.Participation -0.01724555 -0.11284562  0.09265604 -0.02894539
## Importance of components:
##                           PC1    PC2     PC3     PC4     PC5    PC6     PC7
## Standard deviation     2.1360 1.1171 0.87671 0.72027 0.55309 0.5299 0.44935
## Proportion of Variance 0.5703 0.1560 0.09608 0.06485 0.03824 0.0351 0.02524
## Cumulative Proportion  0.5703 0.7263 0.82235 0.88720 0.92544 0.9605 0.98578
##                            PC8
## Standard deviation     0.33728
## Proportion of Variance 0.01422
## Cumulative Proportion  1.00000

As the caption suggests, the data seems to indicate that while participation of women in the labour-force and political decision-making are important components of gender equality, they do not predict the overall (physical) well-being of women. More of this below. The reason why our standardized data is able to provide more insights into this, is due to the standardized means and variations that it contains. PCA creates its principal components (PC) by calculating a line of best fit through the observations of the variables that account for most of the variance (while maintaining the sum-of-squares for the resulting coefficients at =1). The second PC is calculated from the remaining variables similarly, but in such a way that it does not correlate with the first PC. With standardized data, the variations are not affected by the scales of the units of the variables and are as such comparable - no single variable can account for most of the variance due to the large values it takes. Conversely, with the non-standardized data, the variable measuring the per-capita GNI of the countries had the largest explanatory power, since the values that its variation takes are in the tens-of-thousands. The rest of the variables take much lower values (The second highest absolute values are taken by maternal mortality and even that has a range of 0 to 1000), and consequently they seemingly account for less of the variation in the data. This is why we see PC1 correlate almost 100% with the variable GNI.per.Capita in the biplot run on non-standardized variables and why the said variable seems to have a disproportionately high variation among the dataset.

4. Interpreting PC1 and PC2

As indicated above, high labour and political participation of women are not short hand for high level of well being among women. We can see that the data is principally divided between countries where there are high levels death at childbirth and adolescent births, and countries where there are high levels of life-expectancy, education and GNI. We might say that this pricipal component mostly maps reproductive health. Indeed, the former two variables are highly correlated between each other, which is not a surprise: giving birth at a young (teen) age creates a greater risk of death at childbirth due to the fact that the mother’s body is not yet fully developed. The latter four variables on the other hand predict lower levels of maternal mortality and adolescent birth, partly because higher levels of education are known to predict a lower fertility rates, and a delayed start to getting children among women, partly because higher GNI predicts better health care, both contributing to a higher life-expectancy.

On the other hand, labour and political participation presuppose neither of these sets of variables - they are almost completely uncorrelated with them. This indicates that women participate in the labour force and parliament even in poorer countries with more traditional conceptions of child-rearing a worser health-care. On the other hand, it perhaps also indicates that most countries perform suboptimally, when it comes to the participation of women in politics and the labour force. As such those two variables are bad predictors of development either in terms of money, education or health.

5. Multiple Correspondence Analysis on the Tea Data

The tea-dataset, according to the information provided by ?tea, includes three sets of variables: the first 18 relate to tea drinking habits, the next 5 map personal details (with age featuring twice, once as integers and once as a categorical variable with five categories), and the final 13 questions mapping product perception. For this section I have chosen to divide the variables into two sets: the active, habit-variables being fed into the MCA and the supplementary personal detail variables, which will be added as colors into the plot to see if we can find links between these supplementary variables and the two dimensions of variability showcased in the plot. The product-perception variables will be excluded in this analysis. Furthermore, only the categorical age variable will be included.

Below, simple bar graphs have been created alongside the structure and dimensions of the remaining data. The reader is welcome to explor the graphs, but since their interpretation is not a requested part of this analysis, a few short comments should suffice: - We can see that there is a clear distinction between places where tea is drunk (home, friends) and where and when it is not (pub, restaurant, dinner, lunch, work) - People seem to drink Earl Grey frequently - Slightly more women answered the survey.

## [1] 300  36
## 'data.frame':    300 obs. of  36 variables:
##  $ breakfast       : Factor w/ 2 levels "breakfast","Not.breakfast": 1 1 2 2 1 2 1 2 1 1 ...
##  $ tea.time        : Factor w/ 2 levels "Not.tea time",..: 1 1 2 1 1 1 2 2 2 1 ...
##  $ evening         : Factor w/ 2 levels "evening","Not.evening": 2 2 1 2 1 2 2 1 2 1 ...
##  $ lunch           : Factor w/ 2 levels "lunch","Not.lunch": 2 2 2 2 2 2 2 2 2 2 ...
##  $ dinner          : Factor w/ 2 levels "dinner","Not.dinner": 2 2 1 1 2 1 2 2 2 2 ...
##  $ always          : Factor w/ 2 levels "always","Not.always": 2 2 2 2 1 2 2 2 2 2 ...
##  $ home            : Factor w/ 2 levels "home","Not.home": 1 1 1 1 1 1 1 1 1 1 ...
##  $ work            : Factor w/ 2 levels "Not.work","work": 1 1 2 1 1 1 1 1 1 1 ...
##  $ tearoom         : Factor w/ 2 levels "Not.tearoom",..: 1 1 1 1 1 1 1 1 1 2 ...
##  $ friends         : Factor w/ 2 levels "friends","Not.friends": 2 2 1 2 2 2 1 2 2 2 ...
##  $ resto           : Factor w/ 2 levels "Not.resto","resto": 1 1 2 1 1 1 1 1 1 1 ...
##  $ pub             : Factor w/ 2 levels "Not.pub","pub": 1 1 1 1 1 1 1 1 1 1 ...
##  $ Tea             : Factor w/ 3 levels "black","Earl Grey",..: 1 1 2 2 2 2 2 1 2 1 ...
##  $ How             : Factor w/ 4 levels "alone","lemon",..: 1 3 1 1 1 1 1 3 3 1 ...
##  $ sugar           : Factor w/ 2 levels "No.sugar","sugar": 2 1 1 2 1 1 1 1 1 1 ...
##  $ how             : Factor w/ 3 levels "tea bag","tea bag+unpackaged",..: 1 1 1 1 1 1 1 1 2 2 ...
##  $ where           : Factor w/ 3 levels "chain store",..: 1 1 1 1 1 1 1 1 2 2 ...
##  $ price           : Factor w/ 6 levels "p_branded","p_cheap",..: 4 6 6 6 6 3 6 6 5 5 ...
##  $ age             : int  39 45 47 23 48 21 37 36 40 37 ...
##  $ sex             : Factor w/ 2 levels "F","M": 2 1 1 2 2 2 2 1 2 2 ...
##  $ SPC             : Factor w/ 7 levels "employee","middle",..: 2 2 4 6 1 6 5 2 5 5 ...
##  $ Sport           : Factor w/ 2 levels "Not.sportsman",..: 2 2 2 1 2 2 2 2 2 1 ...
##  $ age_Q           : Factor w/ 5 levels "15-24","25-34",..: 3 4 4 1 4 1 3 3 3 3 ...
##  $ frequency       : Factor w/ 4 levels "1/day","1 to 2/week",..: 1 1 3 1 3 1 4 2 3 3 ...
##  $ escape.exoticism: Factor w/ 2 levels "escape-exoticism",..: 2 1 2 1 1 2 2 2 2 2 ...
##  $ spirituality    : Factor w/ 2 levels "Not.spirituality",..: 1 1 1 2 2 1 1 1 1 1 ...
##  $ healthy         : Factor w/ 2 levels "healthy","Not.healthy": 1 1 1 1 2 1 1 1 2 1 ...
##  $ diuretic        : Factor w/ 2 levels "diuretic","Not.diuretic": 2 1 1 2 1 2 2 2 2 1 ...
##  $ friendliness    : Factor w/ 2 levels "friendliness",..: 2 2 1 2 1 2 2 1 2 1 ...
##  $ iron.absorption : Factor w/ 2 levels "iron absorption",..: 2 2 2 2 2 2 2 2 2 2 ...
##  $ feminine        : Factor w/ 2 levels "feminine","Not.feminine": 2 2 2 2 2 2 2 1 2 2 ...
##  $ sophisticated   : Factor w/ 2 levels "Not.sophisticated",..: 1 1 1 2 1 1 1 2 2 1 ...
##  $ slimming        : Factor w/ 2 levels "No.slimming",..: 1 1 1 1 1 1 1 1 1 1 ...
##  $ exciting        : Factor w/ 2 levels "exciting","No.exciting": 2 1 2 2 2 2 2 2 2 2 ...
##  $ relaxing        : Factor w/ 2 levels "No.relaxing",..: 1 1 2 2 2 2 2 2 2 2 ...
##  $ effect.on.health: Factor w/ 2 levels "effect on health",..: 2 2 2 2 2 2 2 2 2 2 ...

Moving on to the actual Multiple Correspondence Analysis. Below, the reader can first find the variable biplot of the MCA, where variables in red denote active variables and variables in green supplementary variables. It is first worth noting that the two dimensions identified by the MCA account together for only 16.7% of the principal inertia, meaning that most of the variation in the data is left unaccounted for. In fact, the summary-statistic shows that increasing the amount of dimensions has little effect on this, with six dimensions only accounting for 36.8% of the principal inertia.

Be that as it may, we will continue to work with these results. The theory that I offer based on the below information is that dimension 1 maps what could be generalised as “tea consumption.” We see the frequency of tea-drinking increase as we move from left to right. Additionally, the variables on the right-hand side of the plot indicate higher levels of tea consumption - to the point that tearooms are visited. The left-hand side showcase answers indicating little interest in tea: the tea is cheap, from the chain store and the people do not generally choose tea as their drink when in restaurants, at home or at friends. Dimension 2 on the other hand maps “tea enjoyment”, or something similar. We see the higher values correspond with upscale tea bought from a teashop in loose-form. The lower variables reflect a more basic approach to tea. But as said, even if correct, these dimensions and their interpretations do not explain a lot of the variance in the data.

Consequently, it is not surprising that when we apply color to the plots to map the respondents by their belonging to different categories in the four supplementary variables, we have a hard time identifying any significant linkages. Age does not seem to be linked to either of the dimensions identified. With sports, there seems to be indication of individuals that do not take part in sport being also lesser tea drinkers, but the indications for this is weak. Similarly, students and elderly seem to drink more tea when it comes to employment status. Perhaps the most clearest linkage is between gender and dimension 1, where females are more prominent at the right-hand side of the scale, perhaps indicating that women tend to be bigger casual tea drinkers, with men more often opting for other drinks, or drinking finer teas, as indicated by the higher amount of men at the top of dimension two.

Overall, these plots are not helpful. Primarily due to their low ability to account for variation in the data, and also, perhaps consequently, due to the absence of other linkages that they might have identified.

## 
## Call:
## MCA(X = Tea_Who_And_How, quanti.sup = NULL, quali.sup = 19:22,  
##      graph = FALSE) 
## 
## 
## Eigenvalues
##                        Dim.1   Dim.2   Dim.3   Dim.4   Dim.5   Dim.6   Dim.7
## Variance               0.148   0.116   0.094   0.078   0.073   0.072   0.067
## % of var.              9.392   7.335   5.972   4.957   4.635   4.534   4.233
## Cumulative % of var.   9.392  16.727  22.699  27.655  32.291  36.824  41.057
##                        Dim.8   Dim.9  Dim.10  Dim.11  Dim.12  Dim.13  Dim.14
## Variance               0.065   0.060   0.059   0.058   0.055   0.052   0.050
## % of var.              4.106   3.829   3.719   3.672   3.500   3.319   3.193
## Cumulative % of var.  45.163  48.992  52.711  56.383  59.883  63.202  66.395
##                       Dim.15  Dim.16  Dim.17  Dim.18  Dim.19  Dim.20  Dim.21
## Variance               0.049   0.046   0.045   0.044   0.039   0.037   0.035
## % of var.              3.094   2.932   2.841   2.785   2.480   2.330   2.210
## Cumulative % of var.  69.488  72.421  75.262  78.047  80.527  82.857  85.067
##                       Dim.22  Dim.23  Dim.24  Dim.25  Dim.26  Dim.27  Dim.28
## Variance               0.034   0.034   0.029   0.028   0.027   0.026   0.024
## % of var.              2.149   2.132   1.860   1.755   1.690   1.637   1.529
## Cumulative % of var.  87.216  89.348  91.208  92.963  94.653  96.290  97.820
##                       Dim.29  Dim.30
## Variance               0.019   0.016
## % of var.              1.196   0.984
## Cumulative % of var.  99.016 100.000
## 
## Individuals (the 10 first)
##                  Dim.1    ctr   cos2    Dim.2    ctr   cos2    Dim.3    ctr
## 1             | -0.573  0.739  0.160 | -0.182  0.095  0.016 | -0.348  0.427
## 2             | -0.381  0.326  0.139 | -0.127  0.047  0.016 | -0.604  1.288
## 3             |  0.095  0.020  0.006 | -0.143  0.059  0.013 |  0.303  0.325
## 4             | -0.634  0.904  0.281 |  0.009  0.000  0.000 |  0.141  0.070
## 5             | -0.138  0.043  0.023 | -0.100  0.029  0.012 | -0.105  0.039
## 6             | -0.734  1.212  0.265 |  0.018  0.001  0.000 | -0.027  0.003
## 7             | -0.084  0.016  0.008 | -0.158  0.072  0.027 | -0.081  0.023
## 8             | -0.262  0.154  0.054 | -0.022  0.001  0.000 | -0.032  0.004
## 9             |  0.198  0.088  0.033 |  0.179  0.092  0.027 | -0.593  1.244
## 10            |  0.344  0.265  0.081 |  0.436  0.546  0.130 | -0.461  0.752
##                 cos2  
## 1              0.059 |
## 2              0.350 |
## 3              0.058 |
## 4              0.014 |
## 5              0.013 |
## 6              0.000 |
## 7              0.007 |
## 8              0.001 |
## 9              0.295 |
## 10             0.146 |
## 
## Categories (the 10 first)
##                   Dim.1     ctr    cos2  v.test     Dim.2     ctr    cos2
## breakfast     |   0.231   0.908   0.049   3.835 |  -0.221   1.062   0.045
## Not.breakfast |  -0.213   0.838   0.049  -3.835 |   0.204   0.980   0.045
## Not.tea time  |  -0.483   3.622   0.181  -7.360 |   0.116   0.269   0.011
## tea time      |   0.375   2.808   0.181   7.360 |  -0.090   0.209   0.011
## evening       |   0.331   1.332   0.057   4.133 |  -0.013   0.003   0.000
## Not.evening   |  -0.173   0.696   0.057  -4.133 |   0.007   0.001   0.000
## lunch         |   0.714   2.656   0.088   5.121 |  -0.352   0.824   0.021
## Not.lunch     |  -0.123   0.457   0.088  -5.121 |   0.060   0.142   0.021
## dinner        |  -0.795   1.568   0.048  -3.769 |   0.827   2.174   0.051
## Not.dinner    |   0.060   0.118   0.048   3.769 |  -0.062   0.164   0.051
##                v.test     Dim.3     ctr    cos2  v.test  
## breakfast      -3.665 |  -0.621  10.329   0.356 -10.315 |
## Not.breakfast   3.665 |   0.573   9.535   0.356  10.315 |
## Not.tea time    1.773 |   0.158   0.607   0.019   2.403 |
## tea time       -1.773 |  -0.122   0.471   0.019  -2.403 |
## evening        -0.169 |   0.510   4.994   0.136   6.383 |
## Not.evening     0.169 |  -0.267   2.611   0.136  -6.383 |
## lunch          -2.520 |   0.356   1.037   0.022   2.552 |
## Not.lunch       2.520 |  -0.061   0.178   0.022  -2.552 |
## dinner          3.922 |   0.637   1.586   0.031   3.023 |
## Not.dinner     -3.922 |  -0.048   0.119   0.031  -3.023 |
## 
## Categorical variables (eta2)
##                 Dim.1 Dim.2 Dim.3  
## breakfast     | 0.049 0.045 0.356 |
## tea.time      | 0.181 0.011 0.019 |
## evening       | 0.057 0.000 0.136 |
## lunch         | 0.088 0.021 0.022 |
## dinner        | 0.048 0.051 0.031 |
## always        | 0.060 0.003 0.053 |
## home          | 0.013 0.000 0.131 |
## work          | 0.099 0.049 0.002 |
## tearoom       | 0.344 0.017 0.000 |
## friends       | 0.208 0.011 0.160 |
## 
## Supplementary categories (the 10 first)
##                  Dim.1   cos2 v.test    Dim.2   cos2 v.test    Dim.3   cos2
## F             |  0.168  0.041  3.511 | -0.114  0.019 -2.375 | -0.045  0.003
## M             | -0.245  0.041 -3.511 |  0.166  0.019  2.375 |  0.066  0.003
## employee      | -0.153  0.006 -1.305 | -0.145  0.005 -1.243 |  0.064  0.001
## middle        | -0.010  0.000 -0.071 |  0.312  0.015  2.114 | -0.330  0.017
## non-worker    |  0.005  0.000  0.045 |  0.175  0.008  1.580 | -0.254  0.017
## other worker  | -0.030  0.000 -0.140 |  0.009  0.000  0.041 |  0.063  0.000
## senior        |  0.370  0.018  2.325 |  0.047  0.000  0.297 | -0.106  0.001
## student       |  0.031  0.000  0.299 | -0.279  0.024 -2.664 |  0.391  0.047
## workman       | -0.454  0.009 -1.601 |  0.216  0.002  0.764 |  0.061  0.000
## Not.sportsman | -0.031  0.001 -0.444 |  0.016  0.000  0.222 | -0.036  0.001
##               v.test  
## F             -0.942 |
## M              0.942 |
## employee       0.545 |
## middle        -2.238 |
## non-worker    -2.285 |
## other worker   0.292 |
## senior        -0.667 |
## student        3.732 |
## workman        0.217 |
## Not.sportsman -0.518 |
## 
## Supplementary categorical variables (eta2)
##                 Dim.1 Dim.2 Dim.3  
## sex           | 0.041 0.019 0.003 |
## SPC           | 0.029 0.044 0.066 |
## Sport         | 0.001 0.000 0.001 |
## age_Q         | 0.008 0.066 0.103 |